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Abstract 

We  derive  new  trace  inequalities  for  NURBS-mapped  domains.  In  addition  to  Sobolev-type 
inequalities,  we  derive  discrete  trace  inequalities  for  use  in  NURBS-based  isogeometric  anal¬ 
ysis.  All  dependencies  on  shape,  size,  polynomial  order,  and  the  NURBS  weighting  function 
are  precisely  specified  in  our  analysis,  and  explicit  values  are  provided  for  all  bounding  con¬ 
stants  appearing  in  our  estimates.  As  hexahedral  finite  elements  are  special  cases  of  NURBS, 
our  results  specialize  to  parametric  hexahedral  finite  elements,  and  our  analysis  also  gener¬ 
alizes  to  T-spline-based  isogeo  metric  analysis. 

Keywords:  trace  inequalities,  isogeometric  analysis,  NURBS,  T-splines,  Nitsche’s  method, 
parametric  hexahedral  finite  elements,  rational  basis  functions 


1  Introduction 

Since  its  introduction  in  2005  by  Hughes,  Cottrell,  and  Bazilevs  [24],  isogeometric  analysis 
has  emerged  as  a  powerful  design-through-analysis  technology.  The  underling  concept  behind 
isogeometric  analysis  is  simple:  utilize  the  same  basis  for  finite  element  analysis  (FEA)  as  is 
used  to  describe  computer  aid  design  (CAD)  geometry.  The  first  instantiations  of  isogeomet¬ 
ric  analysis  were  based  upon  Non-Uniform  Rational  B-Splines  (NURBS),  and  isogeometric 
analysis  has  since  been  extended  to  T-spline  [4,  16]  and  subdivision  [13]  discretizations.  Iso¬ 
geometric  analysis  has  shown  much  promise  in  a  wide  variety  of  application  areas  such  as 
structural  analysis  [15,  17,  28],  fluid-structure  interaction  [5],  and  electromagnetism  [12], 
and  isogeometric  analysis  now  has  a  firm  theoretical  basis  [6,  11,  21],  However,  a  number 
of  research  challenges  remain.  Among  these  challenges  is  the  imposition  of  strong  bound¬ 
ary  conditions.  Unlike  Lagrange  finite  elements,  NLIRBS,  T-spline,  and  subdivision  basis 
functions  do  not  interpolate  function  values  at  nodal  points.  Imposition  of  strong  bound¬ 
ary  conditions  is  even  less  straight-forward  for  fourth-order  systems  (such  as  those  arising 
in  Cahn-Hilliard  phase-held  models  [22]  and  Kirchhoff  plate  and  shell  theories  [25])  where 
strong  boundary  conditions  involving  normal  derivatives  naturally  arise.  Hence,  special  care 


1 


must  be  taken  in  order  to  ensure  an  isogeometric  discretization  satisfies  prescribed  strong 
boundary  conditions. 

Recently,  Nitsche’s  method  [29]  has  been  explored  as  a  means  of  weakly  enforcing  strong 
boundary  conditions  in  isogeometric  analysis  [8,  18].  This  procedure  has  been  found  to 
be  especially  powerful  in  isogeometric  flow  simulation  where  a  close  connection  with  so- 
called  turbulent  wall  models  has  been  discovered  [7].  In  Nistche’s  method,  a  stabilization 
parameter  appears  which  must  be  chosen  large  enough  to  ensure  coercivity.  Notably,  the 
parameter  should  be  chosen  large  enough  such  that  a  discrete  trace  inequality  is  satisfied. 
This  requirement  results  in  Nitsche  parameters  which  necessarily  depend  locally  on  the 
shape  of  the  element  and  the  chosen  discretization,  but  such  dependencies  have  not  yet  been 
made  explicit  in  the  context  of  isogeometric  analysis.  In  this  paper,  new  trace  inequalities 
for  NURBS-based  isogeometric  analysis  are  derived  where  all  dependencies  on  shape,  size, 
polynomial  order,  and  the  NURBS  weighting  function  are  precisely  specified.  Furthermore, 
explicit  values  are  provided  for  bounding  constants  appearing  in  our  estimates  and,  as  such, 
these  inequalities  can  be  directly  utilized  in  the  design  of  Nitsche’s  stabilization  parameter. 

An  outline  of  this  paper  is  as  follows.  In  Section  2,  we  briefly  introduce  requisite  ter¬ 
minology  and  definitions.  In  Section  3,  we  derive  a  new  Sobolev  trace  inequality  for  Hl 
functions  defined  on  NURBS-mapped  domains.  In  Section  4,  we  derive  discrete  trace  in¬ 
equalities  specifically  for  use  in  NURBS-based  isogeometric  analysis.  Finally,  in  Section  5, 
we  draw  conclusions.  Before  proceeding,  we  would  like  to  make  a  couple  of  remarks.  First  of 
all,  the  focus  of  this  paper  will  be  exclusively  the  three-dimensional  setting.  Similar  results 
to  those  presented  here  hold  in  the  two-dimensional  setting,  albeit  with  different  constants 
and  dependencies.  Second,  while  we  only  discuss  NURBS-based  isogeometric  analysis,  our 
results  immediately  apply  to  any  isogeometric  discretization  technique  based  on  rational 
functions.  Notably,  they  apply  to  T-spline  and  parametric  hexahedral  finite  element  dis¬ 
cretizations.  Unfortunately,  our  results  do  not  extend  to  subdivision  discretizations  as  they 
are  not  rational. 


2  Preliminaries 


Throughout  this  paper,  we  make  use  of  the  classical  Lebesgue  spaces  Lg(D )  endowed  with 
norm  ||  •  ||  li(D)  where  1  <  q  <  oo  and  D  C  ld  is  a  generic  open  domain  for  integer  d  >  1. 
We  will  also  utilize  the  Sobolev  spaces  Hk(D )  for  k  a  positive  integer,  endowed  with  norm 
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L2{D). 
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(2) 


We  denote  by  <9 D  the  boundary  of  an  open  domain  D  C  Md,  and  we  further  denote  by 
II  '  ||  L2{dD)  the  Lebesgue  norm  of  order  2  on  <9 D.  We  will  occasionally  employ  the  Sobolev 
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spaces  Wk,00(D)  for  k  a  positive  integer,  and  we  endow  these  spaces  with  norm 


ll/'ll  Wk,oo  .  — 


Qa  1  Qad 

-  - u 

dx^1  ' ' '  d°dXd 
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and  semi- norm 
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Note  that  these  are  not  the  classically  defined  norms  and  semi-norms  for  the  Wk’°°  spaces. 
For  matrix-valued  functions  M  :  D  — >  Mmxn,  we  define  a  spectral  L°°  operator  norm 


||M||l°°(.D),Z 
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xeKm 


(EI=  i  MijXif 
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L°°(D) 
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where  My  and  xt  denote  component-wise  entries  of  M  and  x  respectively  and  crmax  denotes 
maximum  singular  value.  For  third-order  tensor-valued  functions  T  :  D  — >  Kmxnxo;  we 
define  an  analogous  L°°  operator  norm 


|  L°°(D),l 


sup 

x,y£R" 


En  /v^o  v^m  m  \2\  1/2 
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L°°(D) 
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where  Tl]k  denotes  a  component-wise  entry  of  T.  The  norms  dehned  by  (4)  and  (5)  have  been 
specially  chosen  as  to  simplify  the  analysis  presented  in  this  paper.  The  rest  of  this  section 
will  be  devoted  to  a  brief  introduction  of  univariate  and  multivariate  B-spline  basis  functions, 
NURBS  (non-uniform  rational  B-spline)  basis  functions,  and  the  NURBS  geometrical  map 
F.  A  more  complete  introduction  to  NURBS  and  B-splines  may  be  found  in  [30],  and  for  an 
introductory  text  on  NURBS-based  isogeometric  analysis,  see  [14]. 


2.1  Univariate  B-splines 

For  two  positive  integers  p  and  n,  representing  degree  and  dimensionality  respectively,  let  us 
introduce  the  ordered  knot  vector 


“  :  —  {0  —  6)  •  •  •  5  U+P+l  ~  1} 


(6) 


where 


6  —  £2  6:  ■  ■  •  Cn+p+1- 


Given  S  and  p,  univariate  B-spline  basis  functions  are  constructed  recursively  starting  with 
piecewise  constants  (p  —  0): 


1  if  &  <  £  <  6+1 

0  otherwise. 


(7) 
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For  p  —  1,2,3,...,  they  are  defined  by 


sm 


C  Ci  TDP—  l/£\  I  C»+P+l  C 


(8) 


When  £i+p  —  £j  =  0,  ,  -  %  is  taken  to  be  zero,  and  similarly,  when  A +„+i  —  £i+1  =  0, 
&+P+i-£  js  taken  to  be  zero.  B-spline  basis  functions  are  piecewise  polynomials  of  degree 

Si+p+1  Si+1 

p,  form  a  partition  of  unity,  have  local  support,  and  are  non-negative. 

Let  us  now  introduce  the  vector  £  =  {(j, . . . ,  £m}  of  knots  without  repetitions  and  a 
corresponding  vector  {rlf . .  .  ,rm}  of  knot  multiplicities.  That  is,  rl  is  defined  to  be  the 
multiplicity  of  the  knot  Q  in  S.  By  construction,  Y2iLi  ri  —  n  +  p  +  1.  We  assume  that 
rt  <  p  +  1.  Let  us  further  assume  throughout  that  f\  —  rm  —  p  +  1,  i.e,  that  5  is  an 
open  knot  vector.  At  the  point  Q,  B-spline  basis  functions  have  a.j  :=  p  —  r3  continuous 
derivatives.  Therefore,  —  1  <  a3  <  p—  1,  and  the  maximum  multiplicity  allowed,  r3  —  p+  1, 
gives  a  discontinuity  at  ( 3 .  By  construction,  ay  =  am  =  —  1. 


2.2  Multivariate  Tensor-Product  B-splines 

The  definition  of  multivariate  B-splines  follows  easily  through  a  tensor-product  construc¬ 
tion.  Let  us  focus  on  the  three-dimensional  case.  Notably,  let  us  consider  the  unit  cube 
£2  =  (0,  l)3  C  M3,  which  we  refer  to  as  the  patch.  Mimicking  the  one-dimensional  case,  given 
integers  pd  and  rid  f°r  d  =  1,2,  3,  let  us  introduce  open  knot  vectors  . . . ,  £ nd+Pd+id } 

and  the  associated  vectors  Cd  =  {Cm,  •  •  • ,  C md,d},  {n,d,  •  •  • ,  rmd4},  and  ad  =  {altd, amd4}. 
There  is  a  parametric  Cartesian  mesh  Qh  associated  with  these  knot  vectors  partitioning  the 
parametric  domain  £2  into  rectangular  parallelepipeds.  Visually, 

Qh  =  {Q  =  $$d=l,2,3  (C id4i  Cid+l,(i)  ,  1  <  id  <  ™-d  ~  1}  •  (9) 

For  each  element  Q  G  Qh  we  associate  a  parametric  mesh  size  Hq  =  /iQjmax  where  hgm^ 
denotes  the  length  of  the  largest  edge  of  Q.  Also,  for  each  element,  we  define  a  shape 
regularity  constant 


Aq  =  (10) 

min 

where  IiQ>rn i„  denotes  the  length  of  the  smallest  edge  of  Q. 

We  associate  with  each  knot  vector  (d  =  1,2,3)  univariate  B-spline  basis  functions 
Bfj  of  degree  pd  for  i  =  1, . . .  ,  nd.  On  the  mesh  Qd,  we  define  the  tensor-product  B-spline 
basis  functions  as 


Bfif' "  :=  BftBftBft,,  i  =  1, . . . ,  nu  j  =  1, . . . ,  n2,  k  =  1, . . . ,  n3. 


(11) 


Like  their  univariate  counterparts,  multivariate  B-spline  basis  functions  are  piecewise  poly¬ 
nomial,  form  a  partition  of  unity,  have  local  support,  and  are  non-negative. 
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2.3  NURBS 


Let  us  define  a  set  of  positive  weights  Wijk  for  1  <  i  <  n\,  1  <  j  <  n2, 1  <  k  <  ns-  The 
so-called  NLIRBS  (Non-LIniform  Rational  B-Spline)  basis  functions  are  defined  to  be 


TDP1P2P3 


w 


=  l,...,ni,  j  =  l,...,n2,  k  =  1,. . .  ,n3, 


(12) 


where  BP^P2P3  are  the  B-spline  basis  functions  defined  in  the  previous  subsection  and  w  is 
the  weighting  function  defined  by 


w  = 


ni,ri2,n3 

E  srs 


i=l,j‘=l,fe=l 


(13) 


Note  that  as  multivariate  B-spline  basis  functions  form  a  partition  of  unity,  the  NLIRBS 
basis  reduces  to  the  B-spline  basis  when  the  weights  are  chosen  to  be  unity.  Also  note  that 
the  NLIRBS  basis  functions  are  pointwise  positive  and  form  a  partition  of  unity. 

In  NLIRBS-based  isogeometric  analysis,  the  physical  domain  is  defined  through  a  NLIRBS 
geometrical  mapping 

ni,ra2,n3 

f  =  xi  fv-Rgr*  <14) 

i=l,j=l,k=l 

where  P ^  G  M3  are  the  so-called  control  points.  F  is  a  parametrization  of  the  physical 
domain  of  interest,  that  is, 

F  :  (0,1)3  — >■  Cl. 

We  assume  throughout  that  F  is  invertible,  with  smooth  inverse  F^1,  on  each  element 
Q  G  Qh-  NLIRBS  are  capable  of  representing  all  conic  sections,  such  as  circles  and  ellipses, 
and  cylinders,  spheres,  tori,  ellipsoids,  are  also  exactly  representable.  In  general,  more 
than  one  patch  must  be  utilized  to  represent  a  chosen  geometry.  In  this  situation,  proper 
constraints  should  be  enforced  in  order  to  ensure  global  continuity.  For  more  on  multi-patch 
geometries,  see  Chapter  2  of  [14]. 

For  each  element  Q  in  the  parametric  domain  there  is  a  corresponding  physical  element 
K  =  F  (Q).  We  define  the  physical  mesh  to  be 


JCh  =  {K  :  K  =  F(Q),  Q  G  Qh }  =  F (Qh). 
We  define  for  each  physical  element  K  G  JCh  a  physical  mesh  size 

hie  =  ||  VF \\L°°(Q),ihQ 


(15) 

(16) 


where  Q  is  the  pre- image  of  K  and  VF  is  the  matrix  of  partial  derivatives  of  the  coordinate 
components  of  F  (i.e.,  (VF)^.  =  II4).  We  introduce  the  space  Vh  of  NURBS  on  Cl  (which 
is  the  push-forward,  of  the  space  of  NLIRBS  on  the  patch  (0,  l)3) 


Vh  :=  span  {A™3  o  F 


_1  ^  ni,ri2,n3 

f  i= l,j=l,k=l 


(17) 


In  isogeometric  analysis,  the  space  Vh  is  utilized  as  the  trial  space  in  a  Galerkin  or  Petrov- 
Galerkin  method.  Note  that  the  functions  in  Vh  are  smooth  on  each  physical  element  K  G 
JCh,  and  the  NLIRBS  basis  functions  and  weighting  function  are  smooth  on  each  parametric 
element  Q  G  Qh- 
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3  A  Sobolev  Trace  Inequality  for 
NURBS-Mapped  Domains 

Sobolev  trace  inequalities  play  an  important  role  in  the  analysis  of  partial  differential  equa¬ 
tions  and  their  numerical  solution.  Generally  speaking,  these  inequalities  involve  bounding 
the  trace  norm  of  a  function  by  a  constant  multiplied  by  its  interior  norm.  The  precise 
value  of  this  constant  and  its  dependency  on  domain  shape  and  size  has  been  the  focus  of  a 
collection  of  papers  in  the  mathematics  literature.  Optimal  Sobolev  trace  constants  for  the 
half-plane  were  derived  by  Beckner  [10]  and  Escobar  [20]  based  on  earlier  work  of  Lions  [27]. 
Optimal  constants  for  general  domains  were  derived  in  [3],  but  the  derived  constants  are 
specified  in  terms  of  an  awkward  and  computationally  prohibitive  limiting  procedure.  Ex¬ 
plicit,  but  not  necessarily  optimal,  Sobolev  trace  constants  for  simplices  and  parallelepipeds 
were  presented  by  Vesser  and  Verfiirth  in  [32]  for  use  in  finite  element  a  posteriori  error 
estimation.  In  this  section,  we  present  a  new  Sobolev  trace  inequality  for  NURBS-mapped 
domains.  This  trace  inequality  is  completely  explicit  with  respect  to  mesh  size,  parametric 
shape  regularity,  and  local  measures  of  the  NLIRBS  parametric  mapping  F .  and  appropriate 
values  are  provided  for  the  bounding  constant  appearing  in  our  estimate.  It  should  be  noted 
that  while  the  derived  inequality  is  not  necessarily  sharp,  special  care  is  taken  to  ensure  that 
the  estimate  is  easily  computable  for  utilization  in  numerical  simulation. 

3.1  General  Lipschitz  Domains  and  the  Unit  Cube 

We  begin  with  the  standard  Sobolev  trace  theorem  for  general  Lipschitz  domains.  Its  proof 
may  be  found  in  [1],  for  example. 

Theorem  3.1.  If  D  C  M3  is  a  Lipschitz  domain,  then  Hl(D)  is  continuously  embedded  in 
L2(dD).  That  is,  there  exists  a  positive  constant  C*(D)  such  that,  for  all  f  G  H1(D), 

ll/ll!2(aD)  —  C*m\f\\2rn(Dy  (18) 

The  constant  C*(D)  appearing  above  necessarily  depends  on  the  domain  D ,  but  the  the¬ 
orem  does  not  explicitly  specify  such  a  dependence.  In  order  to  arrive  at  such  a  dependence 
for  NLIRBS-mapped  domains,  let  us  first  present  the  following  corollary  to  Theorem  3.1. 

Corollary  3.1.  There  exists  a  positive  constant  Cu  such  that,  for  all  f  G  1L1(( 0,  l)3), 

ll/lli2(<9(o,i)3)  <  Gi || /|| // 1 ((o,i)3)-  (19) 

Note  that  the  optimal  constant  appearing  in  (19)  is  precisely  the  maximizer  of  the 
Rayleigh  quotient 

ll/lll/2(c>(0,l)3)  /  2m 

WfWmmm  [  ’ 

where  /  G  i/1((0,l)3.  Hence,  the  optimal  constant  is  equivalent  to  the  largest  eigenvalue 
A max  °f  the  generalized  variational  eigenproblem:  Find  /  G  L71(( 0,  l)3),  A  G  M+  such  that 

(/>  ^)l2(9(o,i)3)  ~  A  (/,  5,)ifi((o,1)3) ,  V<?  G  H  ((0,1)  ).  (21) 
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In  order  to  arrive  at  the  value  of  \max,  we  have  numerically  solved  (21)  using  Galerkin’s 
method  in  conjunction  with  a  sequence  of  nested  tri-cubic  B-spline  spaces.  These  computa¬ 
tions  reveal  that  the  value  of  Amax,  up  to  ten  significant  figures,  is 

Xmax  =  6.165748028...  (22) 

This  value  provides  a  sharp  lower  bound  for  the  constant  Cu  appearing  in  Corollary  3.1.  It 
is  interesting  that  the  Rayleigh  quotient  (20)  is  nearly  maximized  for  constant  /.  This  is 
also  true  for  the  unit  line  and  the  unit  square. 


3.2  The  Rectangular  Parallelepiped 


We  now  derive  a  Sobolev  trace  inequality  for  the  rectangular  parallelepiped  by  utilizing  a 
scaling  argument.  We  note  that  the  trace  inequality  is  similar  in  nature  to  that  of  Corollary 
4.5  of  [32],  but  our  method  of  proof  and  final  result  are  ultimately  different. 

Lemma  3.1.  Let  D  C  M3  denote  a  rectangular  parallelepiped  whose  length,  width,  and  height 
are  hi,  h2,  and  h3  respectively.  Then,  for  all  f  G  Hl(D), 

II/IIl2(9d)  ^  ' y ~  (ll/lll^)  +  -^l/lre(r>))  (23) 

where  Cu  is  a  positive  constant  chosen  large  enough  such  that  inequality  (19)  of  Corollary 
3.1  is  satisfied  and 

M  :=  max  {h\h2,  h2h3,  hih3}  , 

N  :=  max  {/r2,  /i|,  K\ }  , 

V  :=  hih-2h3  =  the  volume  of  D. 

Proof.  Let  /  G  Hl(D)  denote  an  arbitrary  function.  We  write 

WfWhpD)  =  I  fds. 

JdD 

We  now  expand  the  surface  integral.  We  write 


f  fds  =  [  +  [  +/+/  +/  +/  f2ds  (24) 

JdD  Jd 

top  ”  ^bottom  ^  ^left  ”  bright  ”  ^front  ”  D  back 

where  Dtop  denotes  the  top  of  the  parallelepiped  D  and  the  other  notation  follows.  Before 
proceeding,  let  g  :  [0,  l]3  — »  D  denote  the  unique  face-preserving  affine  function  mapping 
the  (closed)  unit  cube  onto  D,  the  closure  of  the  open  domain  D.  That  is,  g  maps  the  top 
face  of  the  unit  cube  to  the  top  face  of  D  and  so  on.  Let  us  now  parameterize  the  surface 
integral  using  g.  We  write: 


/  fds 

**  -^bottom 


1  rl 


h  i  h2 


h]h2 


h2h3 


' o  Jo 
"1  r  1 


'0  ^0 
"1  r  1 


0  ^0 


/2( g(xi,x2,  l))dx1dx2 


f(g(x1,x2,0))dxidx2 
/2(g(0,  x2,  x3))dx2dx3 
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1  rl 


f2ds  —  h-2h3 


id 


right 


f2(g(l,x2,x3))dx2dx3 


'o  Jo 
-1  rl 


f2ds  =  h\h3 


-^front 


f2(g(x1,0,x3))dx1dx3 


'o  Jo 
"1  rl 


f2ds  =  h\h3 


-^back 


/2(g(xi,  l,x3))dxldx3. 


o  Jo 


Noting  the  double  integrals  in  the  above  expressions  correspond  to  integration  of  the  pullback 
F  —  f  o  g  over  the  faces  of  the  unit  cube,  we  may  insert  the  above  six  equations  into  (24)  to 
obtain 


L2  (<9D)  —  M 


F2dS  =  M||F|||2(a(0)1)3 


where 

By  Corollary  3.1, 


'9(0, l)15 

M  =  max  {h\h2,  h-2h3,  h\h3} 


(25) 


II-^IIl2(9(0,1)3)  <  Cu  (  ||-F|Il2((0,1)3)  +  I^I//1((0,1)3) 


(26) 


where  Cu  is  a  positive  constant  chosen  large  enough  such  that  (19)  is  satisfied.  Applying  a 
change  of  variables  using  the  inverse  of  g  to  map  back  from  the  unit  cube  to  D,  we  obtain 


ll^lli2((o,i)3)  =  ^ll/lli2(D) 


(27) 


and 


h2  h'2  h2  N 

\F\rnmm  =  yWWhiD)  +  y\\dyf\\biD)  +  ^\\dzf\\lHD)  <  (28) 

where  V  =  h\h2h3  and  N  =  max  {h\,  h|,  h2}.  The  proof  follows  from  combining  (25)- 
(28).  ~  □ 

As  a  direct  corollary  of  the  above  lemma,  we  have  the  following. 

Corollary  3.2.  Let  D  C  M3  denote  a  rectangular  parallelepiped  whose  length,  width,  and 
height  are  h\,  h2,  and  h3  respectively.  Let  ho  =  max  {hi,  h2,  h3}  and  let  \d  >  1  denote  the 
local  shape  regularity  constant 

Xd  =  •  1  V  (29) 

nun  {hi,  h2,  h3\ 

Then,  for  all  f  E  Fll(D), 

II/IIl2(9D)  —  Cu\d  (hD  ||/||i2(D)  +  ho\f\H1(D) 


(30) 


where  Cu  is  a  positive  constant  chosen  large  enough  such  that  (19)  is  satisfied. 


3.3  NURBS-Mapped  Domains 

Armed  with  the  theoretical  results  appearing  in  the  previous  subsection,  we  can  finally  prove 
a  trace  inequality  for  H 1  functions  living  on  the  physical  mesh  /C/-(  which  is  explicit  with 
respect  to  size  and  shape.  The  inequality  is  a  simple  consequence  of  pull-back  operations 
between  physical  space  and  parametric  space.  Before  proceeding,  we  would  like  to  remark 
that  the  trace  inequality  provided  here  naturally  extends  to  any  rectangular  parallelepiped 
mapped  using  a  smooth  mapping  F .  Hence,  it  extends  to  iso-  and  subparametric  hexahedral 
finite  elements  which  are  commonly  used  in  practice. 

Theorem  3.2.  Let  K  G  ICh  and  let  Q  =  F_1(/l).  Then,  for  all  f  G  H1( K), 

II / II L2(i9ic)  —  Cu\q\k  (^hfWfWf^K)  +  ^A'l/lrrqic))  (31) 

where  Cu  is  a  positive  constant  chosen  large  enough  such  that  (19)  is  satisfied,  A q  is  the  local 
shape  regularity  constant  of  Q, 

A/c  =  ||cof(VF)||Loo(Qy||det(VF  ^IIl^aoIIVFII^q^,  (32) 

cof(VF)  is  the  cofactor  matrix  o/VF,  and  det(VF_1)  is  the  determinant  o/ V F_1. 
Proof.  Let  /  G  Hl(K).  We  write 


2 

L2(dK) 


fds. 


Let  us  separate  the  surface  integral  into  a  collection  of  integrals  over  “faces” .  That  is,  let  us 
write 


fds  = 


dK 


+ 


Ki 


top 


+ 


-^■bottom 


+ 


-frreft 


+ 


-bright 


-^-front 


fds 


-^-back 


where  we  have  defined,  for  example,  Ktop  =  F  (Qtop)  where  Qtop  is  the  top  face  of  the 
rectangular  parallelepiped  Q.  Note  that,  by  construction, 


Q  <S)rf=l,2,3 (Cid ,rf 5  Cid+l,rf) 


for  some  1  <  id  <  m(i  —  1  where  d  —  1,  2,  3.  We  can  then  write  the  surface  integral  for  the 
top  face  of  K  as 


(F  (£1,62,  Cis+1,3))  <^(£l,£2,  Ci3+l,3)d^ld^2 


where 


d£i  df2 


As  D F  is  smooth  on  Q,  we  have  that 


^(£l>  ^2,  Cis+1,3)  <  l|cof(VF)||Loo(Qy 
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for  all  £1  G  (Cn,i?  Cii+i,i)  and  £2  G  (Ci2,2,  (12+1,2)  where  cof(VF)  is  the  cofactor  matrix  of 

VF.  Thus, 


f'2ds  <  ||cof(VF)||LOo(Q)), 


(F  ((1,  (2,  Cs3+i,i))  d£,id£,2- 


Note  that  the  remaining  double  integral  in  the  above  expression  is  precisely  the  surface 
integral  for  the  square  of  the  pullback  /  o  F  over  the  top  face  of  Q.  If  we  repeat  the  above 
process  for  the  other  faces  of  K  and  sum  all  of  the  resulting  expressions,  we  obtain  the 
inequality 

\\f\\L2(dI<)  ^  ||cof(VF)|U„w),,||/oF||212TO.  (33) 

Now  let  us  apply  the  trace  inequality  given  by  Corollary  3.2.  This  gives 


11/  0  F\\l2{9Q)  <  Cu\q  (V 11/  o  F  |||1(0)  +  hQ\f  o  F  |y(Q))  (34) 


where  A q  is  the  shape  regularity  constant  for  element  Q  and  Cu  is  a  positive  constant  chosen 
large  enough  such  that  (19)  is  satisfied.  An  application  of  change  of  variables,  Holder’s 
inequality  (see  Appendix  A),  and  the  definition  of  the  spectral  norm  for  matrices  results  in 

11/ °  Filing,  <  l|clet(VF-1)||1»m||/||!a(K)  (35) 

l/°  FIh.(O)  <  l|det(VF-,)||I»(K)||VF||i„(<)H|/|5,lm  (36) 

where  det(VF_1)  is  the  determinant  of  VF-1.  Combining  (33)-(36)  gives 

ll/lll»«MO  <  C„A,||oof(VF)|U-(QW  (Cifcyil/lli.m  +  C2fc„|/|ym)  .  (37) 

where 

Ci  =  ||det(VF-1)||LooW 

and 

C2  =  IldetfVF-'JII^IIVFIl!^),,, 
and  combining  (37)  with  the  definition 

h-K  =  II  VF||Loo(Qyi/lQ 

finally  gives  the  desired  result.  □ 


It  should  be  noted  that  the  constant  XqXr  appearing  in  the  above  theorem  is  a  di¬ 
mensionless  measure  of  physical  shape  regularity.  This  can  be  seen  by  observing  that 
||cof(VF)||ioo(Q);/  is  a  measure  of  how  planes  in  Q  are  expanded  under  F.  || VF 
is  a  measure  of  how  lines  in  Q  are  lengthened  under  F,  and  ||det(VF_1)||£,oo(x)  is  an  inverse 
measure  of  how  volumes  in  Q  are  expanded  under  F .  The  multiplication  of  these  three 
quantities  by  the  local  parametric  shape  regularity  gives  the  physical  shape  regularity.  Also, 
recall  that  we  have  a  sharp  lower  bound  estimate  for  the  bounding  constant  Cu  appearing  in 
(31).  Hence,  Theorem  3.2  gives  a  practically  computable  trace  constant.  Finally,  note  that 
if  F  is  taken  to  be  the  identity  mapping,  then  Xr  =  1.  If  we  had  utilized  alternative  matrix 
norms  in  our  analysis  such  as  the  Frobenius  norm,  we  would  have  obtained  Xr  >  1. 


Remark.  Technically,  F  is  a  mapping  of  the  open  set  (0,  l)3  and  consequently  not  a  vector 
held.  That  being  said,  VF  is  a  common  abuse  of  notation  meaning  D F  in  the  sense  of 
calculus  on  manifolds  [26]. 
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4  Discrete  Trace  Inequalites  for  Isogeometric  Analysis 

Discrete  trace  inequalities  belong  to  a  special  class  of  inverse-type  inequalities  which  are 
posed  over  discrete  function  spaces.  Discrete  trace  inequalities  play  an  important  role  in 
the  selection  of  stabilization  parameters  for  Nitsche’s  method  [29],  and  they  also  play  a 
critical  role  in  the  design  of  symmetric  interior  penalty  methods  [2,  34],  Heretofore,  the 
derivation  of  explicit  discrete  trace  inequalities  has  been  largely  limited  to  simplex  elements 
[19,  31,  33].  In  this  section,  we  derive  new  explicit  discrete  trace  inequalities  for  use  in 
NURBS-based  isogeometric  analysis  where  all  dependencies  on  shape,  size,  polynomial  order, 
and  the  NURBS  weighting  function  are  precisely  specified.  As  in  the  previous  section,  special 
care  is  taken  to  ensure  that  the  derived  estimates  are  easily  computable  for  use  in  numerical 
simulation. 


4.1  Discrete  Trace  Inequalities  for  Tensor-Product 
Polynomials  on  Rectangular  Parallelepipeds 


We  begin  this  section  by  deriving  discrete  trace  inequalities  for  tensor-product  polynomials 
on  rectangular  parallelepipeds.  In  order  to  conduct  such  a  derivation,  we  will  need  the 
following  lemma. 

Lemma  4.1.  Let  Vh  =  VP1'P2,P3  denote  the  space  of  tensor-product  polynomials  of  order 
(PuP2,P3)  defined  on  M3.  Then,  for  all  ph  G  Vh , 

11/11^(8(0,1)3)  <Cinv(pllp2lp3)\\ph\\2L2mm  (38) 

where 

Cinv(pi,P2,P3 )  —  2  ((pi  +  l)2  +  {jp2  +  l)2  +  ip 3  +  1)")  •  (39) 

Proof.  Let  ph  G  Vh.  Since  ph  is  a  tensor-product  polynomial  of  order  (pi,p2,P3),  we  can 
write 

PI  ,P2  ,P3 

ph{x1,x2,x3)  =  E  CijkNi(x1)Nj{x2)Nk{x3) 

i=0,j=0,k=0 

where  Nj  is  the  normalized  shifted  Legendre  polynomial  of  order  l  such  that 

l' Nl{x)Nl{x)dX  =  [  l  i{l=j 
and  Cijk  G  M  are  appropriately  chosen  constants.  By  construction, 


P  1,P2  ,P  3 


II/IIl2((o,i)3)  —  ^2  c ijk ■ 


i=0,j=0,k=0 


Now  let  us  decompose 


1  rl  /  P1,P2,P3 


ll/ll22(a(o,i)3)  — 


0  ^0 


y  CijkNi^NjMNkixs)  dx2dx3  + 


i=0,j=0,k=0 
1  rl  /  Pl,P2,P3 


0  ^0 


E  CijkNii^Njix^Nkixz)  dx2dx3+ 


\i=0,j=0,k=0 
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2 


1  /•!  /  P  1,P2,P3 


0  ./O 


y,  C'ijfciVi(xi)iVi(0)iVfc(x3)  )  dxidxz  + 

\i=0,j=0,k=0 


1  /•  1  /  P  1.P2.P3 


0  ./0 


y  )  dx±dx3  + 

\i=0,j=0,k=0 


1  /•!  /  Pl,P2,P3 


0  ./0 


y  CijkNi (xi)Nj(x2)Nk(0)  )  GteiGfc2  + 

\i=0,j=0,k=0 


1  /•!  /  P  1,P2,P3 


0  ./O 


y  CijkNi (xi)Nj(x2)Nk(l)  )  dxldx2. 

\,i=0j=0,fe=0 


(40) 


/•I  />1  /  P1,P2,P3 

'0  Jo 


If  we  take  advantage  of  the  orthonormality  properties  of  our  shifted  Legendre  polynomials, 
we  can  write 

’1>P2,P3  \  2  P2,P3  /  pi  x  2 

y  Cijfc^V'<(0)jy7-(a:2)^Vfc(a?3)  J  dx2dx3  =  y  f  y  CqftW(O) 

^i=0,j=0,/o=0  /  j=0,fe=0  \  0  / 

Now  we  employ  Theorem  3.1  of  [33],  which  states  that 

1/(0)  <(p+  1)11/11^(0,!)  (41) 

for  every  polynomial  /  of  degree  p.  Since  CljkN%  is  a  polynomial  of  degree  p\  for  every 
j,  k,  we  can  utilize  Theorem  3.1  of  [33]  and  orthonormality  to  obtain 

P2,P3  /PI  \  2  P2,P3  pi 

E  E^.(o)  <  (Pl  +  i)2  e  /  e  (CijkNi(xi))2  dxi 

J=0,fc=0  \i=0  /  j=0,fc=0^°  i=0 


i7=U,fc=U 
Pl,P2,P3 


2 

ijk 


=  (Pi  +  l)2  X  cr 

i=0,j=0,A;=0 

=  (Pi  +  i)2llPhllia((o>i)3) 

By  repeating  this  process  to  bound  the  other  five  terms  appearing  in  (40),  we  obtain  the 
desired  expression.  □ 


We  would  like  to  mention  that  while  the  constant  Cmv  appearing  in  the  above  trace 
inequality  is  not  necessarily  sharp,  it  is  completely  explicit  with  respect  to  polynomial  order. 
We  believe  this  to  be  a  significant  advantage  of  our  derived  estimate.  Alternatively,  one  can 
solve  the  generalized  eigenvalue  problem  inferred  by  (38)  to  arrive  at  a  sharper  bounding 
constant  (see,  for  example,  [23]). 


Lemma  4.2.  Let  D  C  M3  denote  a  rectangular  parallelepiped  whose  length,  width,  and  height 
are  hi,  h2,  and  h3  respectively,  and  define  Vh  =  'PPl,P2,P3  to  be  the  space  of  tensor-product 
polynomials  of  degree  (pi,P2,Ps)-  Let  ho  =  max{/ii,  h2,  h3}  and  let  Xp  >  1  denote  the  local 
shape  regularity  constant 


X  d  — 


ho 

min  { hi,  h2,  h3}' 


(42) 
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(43) 


Then,  for  all  ph  E  Vh , 

11/ III*  (3D)  —  Cinv^Dhp\\ph\\2L2(Dy 

where  Cinv  =  Cinv(p\,p2,pfi)  is  the  positive  constant  defined  by  Equation  (39). 

Proof.  Let  ph  E  Vh.  We  utilize  the  same  scaling  argument  as  in  the  proof  of  Lemma  3.1. 
Notably,  letting  g  :  [0,  l]3  — >  D  denote  the  unique  face-preserving  affine  function  mapping 
the  closed  unit  cnbe  onto  D  (the  closure  of  D)  and  defining  Ph  =  ph  o  g,  we  can  write 

\\ph\ \h(dD)  <  M  [  (Phf  dS  =  M||f>'>||l5(£>(0il)S) 

J  9(0, l)3 


where 


M  =  max  {hih.2,  h2h 3,  h\hz}  . 


By  Lemma  4.1, 

11  11 2  n  11  p^ii2 

\\r  IIl2(9(0,1)3)  —  Lot  Ip  IIl2((0,1)3) 

where  Cinv  =  Cinv(pi,p2,p:i)  is  the  positive  constant  deffiied  by  Equation  (39).  A  simple 
change  of  variables  formula  gives 


\\Ph\\hmm 


1 


11/  III2 (D) 


and  the  desired  result  follows  by  concatenating  all  of  the  above  equalities  and  inequalities.  □ 


4.2  Discrete  Trace  Inequalities  for  NURBS-Based  Isogeometric 
Analysis  with  Locally  Constant  Weighting  Function 

In  this  subsection,  we  derive  discrete  trace  inequalities  for  isogeometric  functions  defined 
on  NLIRBS-mapped  domains  where  the  weighting  function  is  locally  constant.  This  setting 
arises  frequently  in  practice  such  as  when  polynomial-based  B-splines  are  employed  instead 
of  more  general  NLIRBS  discretizations.  Furthermore,  this  setting  makes  for  a  much  more 
straight-forward  analysis  which  will  be  expanded  upon  later.  We  begin  with  the  following 
lemma,  which  bounds  the  element  trace  of  an  isogeometric  function  by  its  interior  L 2  norm. 

Lemma  4.3.  Let  K  E  /Q,  and  Q  =  F-1(/l).  Suppose  that  the  NURBS  weighting  function 
w,  defined  by  (13),  is  chosen  such  that  it  is  constant  over  Q.  Then,  for  all  uh  E  Vh ,  defined 
in  (17), 

ll^/l !2(3A')  —  Cinv^Q^Khii  ll'u/l!2(A')  (44) 

where  Cinv  =  Cinv{p\,p2,pfi)  is  the  positive  constant  defined  by  (39),  X q  is  the  local  shape 
regularity  constant  of  Q,  and  A k  is  the  shape  regularity  constant  defined  by  (32). 

Proof.  Let  uh  E  Vh.  To  begin,  we  employ  a  change  of  variables  from  physical  space  to 
parametric  space  and  utilize  a  similar  argument  to  that  used  in  the  proof  of  Theorem  3.2  to 
arrive  at 

l|/|ll2(3A')  —  ||c°f(VF)||Loo(Q)!;||Mh  O  F||iw 
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Noting  that  by  supposition  uh  o  F  is  a  polynomial  over  Q ,  we  invoke  Lemma  4.2  to  obtain 
the  expression 

\\uh\\2L2(dK)  —  C,m^AQ/lQ1||cof(VF)||i/°o(Q)!/||M/l  O  F|| ^2(q) 

where  Cinv  =  Cinv(p\,p2,p:i)  is  the  positive  constant  defined  by  (39),  A q  is  the  local  shape 
regularity  constant  of  Q,  and  fiQ  is  the  length  of  the  largest  side  of  Q.  Finally,  mapping  back 
to  physical  space,  we  have 

Wuh\\l^dK)  <  C'm^Q/iQ1||cof(VF)||Loc(Q)!Z||det(VF-1)||LocW||'u/l|||2w. 

The  lemma  follows  by  recalling  the  definitions  of  A k  and  hj<.  □ 

Before  proceeding,  we  would  like  to  note  the  similarity  in  form  between  the  inequalities 
appearing  in  Theorem  3.2  and  Lemma  4.3.  Notably,  both  contain  a  bounding  constant, 
a  shape  regularity  factor,  and  mesh  scaling  terms.  Furthermore,  the  bounding  constants 
appearing  in  both  inequalities  are  explicitly  given. 

We  now  present  an  inequality  which  bounds  the  element  boundary  normal  derivative 
of  an  isogeometric  function  by  its  interior  H 1  semi-norm.  Such  an  inequality  plays  an 
important  role  in  the  design  and  analysis  of  Nitsche’s  method  as  applied  to  second-order 
elliptic  problems. 

Lemma  4.4.  Let  K  e  /Q,  and  Q  =  F_1(iL).  Suppose  that  the  NURBS  weighting  function 
w  is  chosen  such  that  is  constant  over  Q.  Then,  for  all  uh  G  Vh , 

\\Vuh  ■  n\\2L2{dK)  <  CBCinv\Q\Khfl\uh\2Hl(K)  (45) 

where  n  is  the  unit  outward-facing  normal,  Cinv  =  Cinv(pi,p2,P3 )  is  the  positive  constant 
defined  by  (39),  A q  is  the  local  shape  regularity  constant  of  Q,  A k  is  the  shape  regularity 
constant  defined  by  (32), 

Cb  =  +  V^H VF ||Loo(q)!;|| V2F  1\\L°c(K),ihK^J  ,  (46) 


and  VJF  1  is  the  third-order  tensor-valued  function 

d2Fr 1 


(V2F-‘),rt  = 


vk  dxjdxk 


(47) 


Proof.  Let  uh  G  Vh.  To  begin,  we  employ  a  change  of  variables  from  physical  space  to 
parametric  space: 

l|V^-n||iW)  <  /  |V^|2ds 

JdK 

<  ||cof(VF)||Loc(W  f  |  (VF_1  o  F)r  V  (uh  o  F )  |2dt.  (48) 

JdQ 

Note  that  we  cannot  immediately  utilize  Lemma  4.2  as  the  components  of  the  vector 

(VF-1  o  F)T  V  (uh  o  F) 
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are  not  necessarily  polynomial.  For  this  reason,  we  employ  the  decomposition 


VF"1  oF=D  =  D  +  D/ 


(49) 


where 

d  =  mc) 

and  £c  is  the  centroid  of  Q.  Inserting  the  decomposition  (49)  into  (48),  applying  the  Cauchy- 
Schwarz  inequality,  invoking  the  definition  of  the  spectral  norm  for  matrices,  and  then  ap¬ 
plying  Holder’s  inequality  on  the  resulting  term  involving  D;,  we  obtain 


\^Uh  ■  n\\L2(dK)  ^  Ccof 


(^Jqq\ d'v^oF)  |2dt)  7  +  \\V,\\Lo0{Q)}l\uhoF\m{dQ)'\ 


where 


CCOf  -  ||C0f(VF)||Loo(Q)jZ 
We  can  now  finally  invoke  Lemma  4.2,  giving  us  the  result 

||Vuh-n|| lHaK)  <  CtraceCcof  Uj^\ D'v^oF)!2^  7  +||D'| 


L°°(Q),l\uh  °  F|ffi(Q) 


(50) 


where 


C' trace  CinviP  \  •  P'2‘  • 

Again  invoking  the  decomposition  (49),  the  Cauchy-Schwarz  inequality,  the  definition  of  the 
spectral  norm  for  matrices,  and  Holder’s  inequality,  we  obtain 

Qf|  D'vta'oFlpisj  <  (  J  |DtV  (uh  o  F)  |2d$|  +  ||D'|U»(W|«'‘  °  F|„iw). 


(51) 


By  concatenating  (50)  and  (51),  we  arrive  at 

\\Vuh-n\\2L2{dK)  <  CtraceCcof  ^ (^Jq  |DTV  (wh  o  F )  \2d(^j  '  +  2||D/| 


L™(Ct),l\uh  0  F|h!(9Q) 


(52) 


Note  that  a  change  of  variables  gives 


T-r t  (~.h 


ID  V 


'Q 


\  1/2 

(uh  o  F)  \2dS)  <  lldet  (VF-1)  W^{K)\uh\^{K) 


(53) 


and  (36)  gives 


u 


°  Flr/pQ)  <  || det  (VF  x)  \\Ll(K)\\'VF\\L°o(Q)7i\uh\Hi(Ky 


(54) 
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To  complete  the  proof,  we  recognize  that,  by  Taylor’s  theorem, 


ID'I 


L°°(Q),l 


D-D 


L°°(Q),l 


\/3 

V3 

<  -^hQ\\VF\\Loo(Q)j\ 


v2f_1i 


L°°(K),l 


(55) 


Combining  (52)-(55)  with  the  definitions  of  the  physical  mesh  size  hx  and  the  physical  shape 
regularity  constant  A k  results  in  the  desired  final  expression.  □ 


Note  that  if  we  take  the  physical  mesh  size  hx  — >  0  and  keep  the  parametric  mapping  F 
fixed,  the  constant  Cb  appearing  in  the  above  theorem  will  tend  to  1.0.  Hence,  we  recover 
a  similar  inequality  to  that  appearing  in  Lemma  4.3.  This  suggests  that  we  may  be  able  to 
ignore  in  practice  higher-order  terms  due  to  the  nonlinear  geometrical  mapping. 

By  employing  a  similar  method  of  proof  to  that  of  Lemma  4.4,  we  can  obtain  discrete 
trace  inequalities  for  higher-order  boundary  derivatives.  For  example,  we  can  show  that  for 
a  given  element  K  G  K-h  with  Q  =  F_1(iF), 

\\^uh\\\2{dK)  ^  CBCinv\Q\xhx\\uh\\2H2(K ) 

for  every  uh  G  Vh  where 

Cb  =  1  +  0{hx)- 

Such  an  inequality  plays  an  important  role  in  the  design  and  analysis  of  Nitsche’s  method 
as  applied  to  fourth-order  elliptic  problems.  It  can  be  shown  that  the  explicit  form  for  the 
constant  Cb  appearing  in  the  above  inequality  is  substantially  more  complicated  than  the 
corresponding  constant  appearing  in  Lemma  4.4,  and  this  trend  only  worsens  with  further 
differentiation.  Moreover,  note  that  a  full  H2  norm  appears  rather  than  just  the  H2  semi¬ 
norm.  However,  in  practice,  it  may  be  sufficient  to  replace  Cb  with  one  and  the  H 2  norm 
with  the  H2  semi-norm. 


4.3  Discrete  Trace  Inequalities  for  NURBS-Based  Isogeometric 
Analysis  with  General  Weighting  Function 


Finally,  we  are  ready  to  derive  discrete  trace  inequalities  for  NLIRBS-based  isogeometric 
analysis  with  general  weighting  function.  We  begin  with  the  following  theorem,  which  is  a 
straight-forward  extension  of  Lemma  4.3  to  the  general  setting. 

Theorem  4.1.  Let  K  G  K,h  and  Q  =  F_1(/l).  Then,  for  all  uh  G  Vh, 

llM  IIl2(9A')  —  CwCinv\Q\KhK  ||rt  \\lHK)  (56) 

where  CJinv  =  Cinv{p\,p2,pf)  is  the  positive  constant  defined  by  (39),  A q  is  the  local  shape 
regularity  constant  of  Q,  \x  is  the  shape  regularity  constant  defined  by  (32),  and 

2 

II'u;IIl°o(q)-  (57) 


w 
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Proof.  The  proof  proceeds  in  a  similar  manner  to  that  of  Lemma  4.3.  Let  uh  G  Vh.  We 
employ  a  change  of  variables  from  physical  space  to  parametric  space  to  arrive  at 


\Wh\\h{dK)  <  ||cof(VF)||  L°°(Q),l\\uh  °  F \\2L2 (QQy 
Now,  we  utilize  Holder’s  inequality  to  obtain 


U  II L2{dK)  <  ||cof(VF)||L<x>(Qy 


w(uh  oF)\\2L2{dQ). 


Noting  that  w[uh  o  F  )  is  a  polynomial  over  Q ,  we  invoke  Lemma  4.2  to  obtain  the  expression 

1 2 


«  IIl2(sk)  <  Cinv\QhQ  ||cof(VF 


L°°(Q) 


w(uho  F)|||2(q) 


where  Cinv  =  Cmv(p\,p2,p-i)  is  the  positive  constant  defined  by  (39),  A q  is  the  local  shape 
regularity  constant  of  Q,  and  hq  is  the  length  of  the  largest  side  of  Q.  Finally,  using  Holder’s 
inequality  again  and  then  mapping  back  to  physical  space,  we  have 


uh\\l 2{dI<)  <  Cw Cinv Aq hg1 1 1 cof ( V F ) 1 1 £oo (q) ^ 1 1 det ( V F  x) |U-(iv),z || || |2 {K) 


where 


cw  = 


w 


L°°(Q) 


L°°(Q) 


The  lemma  follows  by  recalling  the  definitions  of  A  a-  and  hj<- 


We  now  present  an  extension  of  Lemma  4.4  to  the  general  setting. 
Theorem  4.2.  Let  K  G  /Q,  and  Q  =  F-1(/l).  Then,  for  all  uh  G  Vh, 

||Vw  -  n||x,2 (dK)  —  CNCinv^Q^I<hK  \u  \H i(K) 


□ 


(58) 


where  Cinv  =  Cmv(pi,p2,p-i)  is  the  positive  constant  defined  by  (39),  A q  is  the  local  shape 
regularity  constant  of  Q,  A k  is  the  shape  regularity  constant  defined  by  (32), 


cN=(cl/2  +  cl/2)\ 


and 


Cl  =  c 


B 


\L°°(Q) 

c2  =  fvr-1C'(ied|VF-1| 


IHI L°°(Q)  "k  it  ^det\  VF  Wl°°(K),1  I I rvi ,o° (Q)  )  i 


L°°(K),l 


w 


L°°(Q)  h'K  I  ) 


C'B  —  (l  +  V/3||VF||/.™(Q)I;||V2F  1||L™(ft-),;/lA')  ) 

Cdet  =  ||det(VF)||^w)||det(VF-1)||ti(ST 


(59) 

(60) 

(61) 

(62) 

(63) 


17 


Proof.  Let  uh  G  Vh.  We  employ  a  change  of  variables  from  physical  space  to  parametric 
space  to  arrive  at 


\Vuh  ■  n\\2L2{aK)  <  I  \Vuh\2ds 

JdK 

<  ||cof(VF)||io »mi  f  \BTV(uhoF)  \2dt 

JdQ 

<  ||cof(VF)||x,oo(Q))/  f  DTV  (—w[uh  o  F 

JdQ  \w 


where  we  have  denoted 


The  product  rule  gives 


D  =  VF'1  o  F. 


V  ( —w(uh  oF)  |  =  —  V  (w(uh  oFj)  +w(/oF)V  (  - 
\w  J  w  v  '  \w 


which  in  conjunction  with  (64)  and  the  Cauchy-Schwarz  inequality  provides  the  expression 

II  •  n|| \i{qk)  <  (A1/2  +  B1/2)2  (65) 

where 

2 

A=  ||cof(VF)||L=c(Q)),  [  — DTV  (w(uh  o  F))  dt 

JdQ  w 

and 

B  =  ||cof(VF)||Loo(Q);Z  [  w(uh  o  F)DrV  f-)  dt. 

JdQ  \w/ 

Invoking  Holder’s  inequality  and  the  definition  of  the  spectral  norm  for  matrices,  we  can 
write 

||cof(VF)||Loo(W  -  /  |DtV  (w(uh  o  F  ))  |2  dt  (66) 

W  L°°(Q)  J dQ 

and 

1  2 

B<  ||cof(VF)||i»(0),,  -  l|VF-1||l„mi,||TO(«*oF)||J5TO)  (67) 

W  W1’°°(Q) 

where  we  have  taken  advantage  of  our  special  definition  for  the  W1,oc  semi-norm.  We  now 
proceed  on  separate  paths  to  bound  the  A  and  B  terms.  We  begin  with  the  A  term  (66). 
Using  an  argument  identical  to  the  one  appearing  in  the  proof  of  Lemma  4.4,  we  obtain  the 
bound 

1  2 

A  <  CsCinyXQXKhjJ  —  (w°F  1)uI>\2h1(k) 

W  L°°(Q) 

where  Cinv  =  Cmv(p],p2,p:i)  is  the  positive  constant  defined  by  (39),  Xq  is  the  local  shape 
regularity  constant  of  Q,  X k  is  the  shape  regularity  constant  defined  by  (32),  and 

Cb  =  ^1  +  \/3 1|  VF  ||z-°°(Q),z ||  V2F  1\\L°°(K),ihK j  • 
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The  product-rule,  Cauchy-Schwarz,  and  Holder’s  then  give 


-i 


A  <  C'B  Cinv  XqXk  hK 


1 

w 


2 

L°°(Q) 


(68) 


where 

^1  =  \\w\\l°°(Q)  \U  I H1^) 


and 


A-2 


W\w1’°°(Q ) 


IVF 


-H 


(K),l 


U 


1!  w 


We  now  proceed  to  bounding  the  B  term.  Recalling  (67)  and  employing  Lemma  4.2,  we 
have 


B  <  Cinv\QhQ  ||cof(VF)||Loc(Q)iZ 


’  I 

Invoking  Holder’s  inequality,  we  further  write 

1  2 

B  ^  ||  cof(^7  F )  ||  ~ 

W  W^°°{Q) 

We  finally  map  back  to  physical  space,  obtaining  the  expression 

1 2 


\\l°°(k),i\\w(u  °  F)||L2(q). 


w\\l°°{Q)  II^F  \\l°°{k),i\\u  °  F|Il2((5)- 


-i 


B  Cinv  ^Q^i<hK 


w 


L°°(Q) 


1 II  L°°(A'),Z  11^^  II  L2(K)  • 


(69) 


I W1’00^) 

Collecting  Equations  (65),  (68),  and  (69)  and  rearranging  terms,  we  have 

2 


Tuh  ■  n\\2LHdK)  <  Cm,,,\rj\Kh^  (A1'2  +  B'/2)' 


(70) 


where 


A  =  C, 


B 


L°°(Q) 


(|MIl°°(q)  W'Ihhk)  +  \w\wi, oc(Q)  ||VF  1\\L°°(K),l\\uh\\L2(K)Sj 


and 

IMIl~(q)  l|VF_1||Z/oo(A')iZ||ufe||L2(K) 

Our  objective  is  to  bound  the  H-^  ||z,2(A')  terms  in  A  and  B  with  terms  proportional  to 
I uh \h 1(k)-  We  recognize  that 


II  '  nll InaK)  =  II ^  K  -  zh)  ■  nll! *{dK) 

for  all  constant  zh  G  Vh.  Let  us  choose  zh  to  be  equal  to 


z 


h 


(uhoF)d£ 


(71) 
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where  \Q\  is  the  volume  of  the  parallelepiped  Q.  Using  change  of  variables  and  Holder’s 
inequality,  we  have 

Ilf1  -  A\lhk)  <  ||det(VF)||^w)||U'‘  o  F  -  ;'>||t,(0). 

Since  zh  is  equal  to  the  average  of  uh  o  F  over  Q,  we  can  employ  Poincare’s  inequality  for  a 
rectangular  parallelepiped  domain  [9]  to  further  write 

\\uh  -  zh\\L2{K)  <  7r-1hQ||det(VF)||J//i(Q)|u,lo  F|^i(q). 

Finally,  we  map  back  to  K: 

ll«'*  -  A\l HK)  <  x-,ftK||det(VF)||^w)||det(VF-1)||i2(J0|«'>|„,(J0.  (72) 

By  replacing  uh  with  uh  —  zh  in  (70)  and  applying  inequality  (72),  we  obtain  the  desired 
final  expression.  □ 

Note  that  if  we  take  the  physical  mesh  size  hj<  — >  0  and  keep  the  parametric  mapping 
F  fixed,  the  constant  Cn  appearing  in  the  above  theorem  will  tend  to  1 1 iy 1 1 IIwIIl°°(Q)- 
Hence,  we  recover  a  similar  inequality  to  that  appearing  in  Theorem  4.1  with  mesh  refine¬ 
ment.  This  suggests  that  we  may  be  able  to  ignore  in  practice  higher-order  terms  due  to  the 
nonlinear  geometrical  mapping.  Moreover,  if  the  NURBS  weighting  function  is  chosen  such 
that  it  is  locally  constant,  we  recover  the  inequality  appearing  in  Lemma  4.4.  Finally,  by 
employing  a  similar  method  of  proof  to  that  of  Theorem  4.2,  we  can  obtain  discrete  trace 
inequalities  for  higher-order  boundary  derivatives.  For  example,  we  can  show  that  for  a  given 
element  K  e  ICh  with  Q  =  F-1(iF), 


\\^uh\\2L2(dK) 


-1 


A  f  '  .\('i  n  r  ^  1\  ^ 


U 


1 HHK) 


for  every  uh  G  Vh  where 

2 

IIu’IIl”(q)  +  0(Iik)- 

L°°(Q) 


5  Conclusions 

In  this  paper,  we  have  derived  new  trace  inequalities  for  use  in  NLIRBS-based  isogeometric 
analysis.  All  dependencies  on  shape,  size,  polynomial  order,  and  the  NLIRBS  weighting  func¬ 
tion  are  precisely  specified  in  our  analysis,  and  explicit  values  are  provided  for  all  bounding 
constants  appearing  in  our  estimates.  Consequently,  these  inequalities  can  be  directly  uti¬ 
lized  in  the  design  of  stabilization  parameters  appearing  in  Nitsche’s  method,  an  attractive 
candidate  for  the  weak  enforcement  of  strong  boundary  conditions  in  isogeometric  analysis. 
As  hexahedral  finite  elements  are  special  cases  of  NURBS,  our  results  specialize  to  paramet¬ 
ric  hexahedral  finite  elements.  Moreover,  as  our  results  are  local  in  the  sense  that  they  apply 
element-wise,  our  analysis  generalizes  to  T- spline-based  isogeometric  analysis. 
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A  An  Alternate  Holder  Inequality 

Lemma  A.l.  Let  D  C  Rd  denote  an  open  domain  for  d  a  positive  integer.  If  f  G  L°°(D ) 

and  g  G  L2(D),  then 

\\fg\\L2(D)  <  II  f  \  \  L°°  (D)  ||fi,||l/2(Z))  •  (73) 

Proof.  Let  /  G  L°°(D )  and  g  G  L2(D).  By  construction, 

l|/fi'IU2(D)  =  ||/2^Jli2(£))- 

By  the  classical  Holder  Inequality, 

\\fg\\L^{D)  <  (II  /2  ||l°°(D)  Ilf?2  ||  L1(D))  / 

/  \  1/2 
=  (ll/||L°°(i?)ll5,lli2(D)J 

=  ||/||l<-(D)||5,||l2(d)-  (74) 

□ 


23 


